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AIRFOIL PROFILES FOR MINIMUM PRESSURE DRAG AT SUPERSONIC 
VELOCITIES - GENERAL ANALYSIS WITH APPLICATION TO 
LINEARIZED SUPERSONIC FLOW 
By Dean R. Chapman 
SUMMARY 

A theoretical investigation is made of the airfoil profile for mini- 
mum pressure drag at zero lift in supersonic flow. In the first part of 
the report a general method is developed for calculating the profile 
having the least pressure drag for a given auxiliary condition, such as 
a given structural- requirement or a given thickness ratio. The various ‘ 
structural requirements considered include bending strength, bending 
stiffness, torsional strength, and torsional stiffness. No assumption 
is made regarding the trailing-edge thickness; the optimum value is deter- 
mined in the calculations as a function of the base pressure. 

To illustrate the general method, tlie optimum airfoil, defined as 
the air foil, having minimum pressure drag for a given auxiliary condition, 
is calculated in a second part of the report using the equations of 
linearized supersonic flow.' It is found that the optimum airfoil in most 
cases has a blunt trailing edge. It also is found that the optimum 
thickness distribution depends only on one dimensionless parameter, termed 
the "base pressure parameter." This parameter involves the Mach number, 
airfoil thickness ratio, and base pressure coefficient. The effect of 
variations in each of these latter three quantities on the, shape of the 
optimum profile is discussed, and a simple criterion formulated for deter- 
mining the' condition under which the optimum trailing— edge thickness is 
greater than zero. The calculated pressure drag of the optimum profile 
is compared to that of a biconvex sharp-trailing— edge profile satisfying 
the same structural requirement. The reduction in pressure drag depends 
on the base pressure parameter and varies from a few percent to as much 
as 75 percent. A rapid graphical method of determining an optimum profile 
is developed and applied to several illustrative examples. 

INTRODUCTION 

In supersonic flow the finite thickness of an airfoil invariably 
introduces a certain amount of pressure drag which can be minimized by 
a rational choice of airfoil shape. The profile for minimum pressure 
drag depends, among other things, on the particular auxiliary condition 
that is imposed on the airfoil geometry. For example, if it is required 
that the optimum profile ‘ (defined herein as the profile of least pressure 
drag for a given auxiliary condition) satisfy the auxiliary condition of 
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a given thickness ratio, then according to a well-known result of Ackeret's 
linearized airfoil theory, the so-called double-wedge profile represents 
the optimum sharp— trailing-edge airfoil. This particular auxiliary con- 
dition, however, does not represent practical cases where an airfoil must 
satisfy a certain structural requirement, such as a given area moment of 
inertia, or a given section modulus. Consequently, Puckett 1 in this 
country, and later Drougge in Sweden (reference l), have made more' elabo- 
rate theoretical analyses to determine the optimum profile for the auxil*- 
iary conditions of a given bending stiffness of the airfoil, and also for 
the condition of a given torsional' stiffness. Both Puckett and Drougge 
used linearized airfoil theory and considered only sharp^trailing-edge 
airfoils. Their results are somewhat limited in two respects: they do 

not cover cases outside the scope of linearized airfoil theory, and, 
though they include the auxiliary conditions of given bending and tor- 
sional stiffness, they do not include the auxiliary condition of a given 
bending strength (given section modulus). 2 A far more important limita- 
tion of these analyses, though, is the tacit assumption that the optimum 
airfoil will have a sharp trailing edge. 

There is a small amount of experimental evidence in the measurements 
of Ferri (reference 2) on airfoils with sharp trailing edges which sug- 
gests that the optimum airfoil might, in fact, have a moderately thick 
trailing edge. The measured profile drag of one airfoil tested by Ferri 
(G.U. 3 airfoil at a Mach number of 1.85) was considerably lower than 
inviscid theory would indicate. Allowance for skin friction would cause 
this discrepancy to become even greater. Schlieren photographs and 
pressure-distribution measurements showed that viscous effects effectively 
thickened the airfoil shape near the trailing edge. From these results 
it can be inferred that at moderate supersonic velocities it is possible 
for an -airfoil with a thickened trailing edge — that is, a blunt— trailing- 
edge airfoil — to have lower drag than a corresponding sharp— trailing-edge 
airfoil. Employing a different approach, this inference has been obtained 
from quantitative considerations in reference 3, where a reasonable esti- 
mate of the base pressure was made and the drag calculated as a function 
of trailing-edge thickness. 3 Such calculations, though very approximate, 
have indicated that in certain cases a moderate increase in trailing-edge 
thickness will decrease the over-all pressure drag. 

1 Puckett , s investigation, made in 1945, indicated precisely the same 
results as Drougge* s analysis, but apparently has not yet been pub— 
lished in a generally available form. 

‘In Puckett's analysis an approximate solution for this case was formu- 
lated on the assumption that the airfoil surface everywhere had con- 
tinuous curvature. As is shown later, however, this assumption is not 
adequate since the surface of the optimum airfoil for a given section 
modulus has two discontinuities in curvature. 

3 Analogous calculations without consideration of the structural charac- 
teristics recently have been made in England by E. S. Eae. 
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Apart' from the reasons just cited for expecting that the optimum 
supersonic airfoil might have a thick trailing edge, there are other 
independent considerations which suggest the same result. By analyzing 
conditions at infinite Mach number, Saenger pointed out in 1933 that even 
with a vacuum at the base the optimum airfoil for a given thickness ratio 
would, in this extreme case, have a trailing— edge thic kne ss equal to the 
maximum airfoil thickness. (See reference 4.) In reference 5, Ivey 
obtained a similar result by calculating the pressure drag at a Mach 
number of 8 for a family of airfoils having various positions of maximum 
thickness. More recently, Smelt (reference 6) has developed an approxi- 
mate condition determining when an airfoil with maximum thickness at the 
trailing edge has lower drag in hypersonic flow than an airfoil with a 
sharp trailing edge. Saenger, Ivey, and Smelt, however, did not consider 
airfoils having a trailing— edge thickness less than the mavimirm airfoil 
thickness, and hence their results do not determine the optimum profile 
for hypersonic velocities. Nevertheless, it is evident that at high 
supersonic Mach numbers the optimum profile has a relatively thick trail- 
ing edge. On this basis it is not unreasonable to expect that at lower 
supersonic Mach numbers the optimum profile would have some thickness at 
the trailing edge. 

The physical reason why it is possible for a blunt-trailing-edge 
airfoil in supersonic flow to have a lower pressure drag than a corre — 
sponding sharp-trailing-edge airfoil is quite simple, as can be illus- 
trated by the two profiles shown in figure 1. These profiles have the 
same area, which corresponds to the same, torsional stiffness of a thin- 
skin structure. The blunt-trailing-^dge airfoil has a smaller thickness 
ratio and a position of maximum thickness which is farther rearward, 
hence the leading-edge angle is smaller and the pressure drag of the sur- 
face forward of the trailing edge is less than that of the sharp-trailing- 
edge airfoil. A certain amount of base drag, however, obviously is added 
by employing a thick trailing edge. If the added base drag is less than 
the reduction in pressure foredrag, then the net result is a smaller 
total pressure drag for the blunt— trailing-^edge airfoil. This invariably 
is the case at extremely high supersonic Mach numbers where the base- drag 
is negligible compared to the pressure foredrag. At low supersonic Mach 
numbers, though, the base drag can be many times the pressure foredrag, 
and the optimum trailing— edge thickness must be expected a priori to 
depend to a great extent on the base pressure. 

The present theoretical analysis was initiated in view of the fore- 
going considerations. The primary purpose of the investigation is to 
develop a method of determining the supersonic airfoil profile for mini- 
mum pressure drag at zero lift, without making an arbitrary assumption 
about the trailing-edge thickness. The profile so determined, which is 
termed an optimum profile, is considered to depend on the base pressure, 
Mach number, and the particular auxiliary condition imposed on the air- 
foil. A .secondary purpose of the investigation is to develop a method 
of sufficient generality to enable second— order and shock— expansion 
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theories to be used in calculating optimum profiles. Such generality is 
desirable in order to obtain results that- are valid at hypersonic Mach 
numbers. 


NOTATION 


B base, pressure parameter for linearized supersonic flow 

-Pv 


r- 7 Me, 2 ' - 1 


L(t/c) 

B limiting value of the base pressure parameter below which the 

optimum airfoil has a blunt trailing edge 

c airfoil chord 

c^ section pressure drag coefficient 

f symbol for the function Py’ 

F(k,cp) incomplete elliptic integral of the first kind of modulus k and 
amplitude cp 

h trailing-edge thickness 

H dimensionless trailing— edge thickness (h/t) 


given value ’of the auxiliary integral 


y* 


dx 


Z 

L 

M 

n 

P 

P 


_ c Jo (t/2) a . 

length of chord over which airfoil thickness is constant 


dimensionless length of chord over which airfoil thickness is 
constant (z/s) 

Mach, number 

arbitrary parameter appearing in the definition of the auxiliary 
integral I 

(For the examples considered n is taken as 1, 2 > 3 or ».) 
static pressure on airfoil surface 


pressure coefficient 


(]>-pj/ ( ) 


base pressure coefficient 
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s 

t 

V 

x 

X 

y 

T 

7 

a 

p 


o. 

l 

00 

b 

V 


t 


distance from leading edge to first downstream position of 
maximum thickness 

mft-srimum thickness of airfoil 

velocity. 

distance from leading edge 

dimensionless distance from leading edge (x/s) 
ordinate of upper half of airfoil 
dimensionless ordinate [y/(t/2)] 
ratio of specific heats (l.4 for air) 

arbitrary parameter appearing in the definition of the auxiliary 
integral I 

(For the examples considered a is taken as 0 or 1.) 
mass density 

Subscripts 

airfoil surface at leading edge 
airfoil surface at trailing edge 
free stream 
base 

vacuum at base 

Superscript 

differentiation with respect to x 

THEORETICAL ANALYSIS 


Assumptions- and Statement of Problem 

/ 

In the analysis which follows several simplifying assumptions are made 
Two-dimensional airfoils in a purely supersonic flow at' zero lift only are 
considered. It is assumed that the pressure at any point on the airfoil 
. surface forward of the trailing edge can be calculated from the flow of an 
inviscid, nonconducting gas. It is further assumed that the leading edge 
is sharp. No analogous assumption is made regarding the trailing edge, 
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however, and consequently a parameter determining the base drag enters 
the problem. Thus, instead of assuming that the- trailing edge is sharp, 
as is usually done, it is assumed that the base pressure coefficient P^, 
is known. This. enables the optimum trailing— edge thickness to be calcu- 
lated as a function of P^, but means that experimental data on base 
pressure in two-dimensional flow are required in order to apply the theo- 
retical results of the analysis. 

Subject to the foregoing assumptions, the problem formulated is that 
of calculating the optimum profile which has the least possible pressure 
drag under a given auxiliary condition, such as a prescribed structural 
requirement, or a given thickness ratio. From the fact that the surface 
pressures on the top and bottom of an airfoil can be calculated independ- 
ently in a supersonic flow, it follows that at zero lift the optimum 
profile will be symmetrical about the chord plane. Consequently, refer- 
ence is made throughout to the thickness distribution of only the upper 
surface of an optimum profile. 

In comparing the pressure drag of various profiles, the chord length 
is held constant, and the thickness distribution along the chord is varied 
in a manner which is arbitrary except for the requirement of satisfying 
the particular auxiliary condition being considered. The various auxil- 
iary conditions investigated are: a given torsional stiffness of the air- 
foil section, a given torsional strength, a given bending stiffness, a 
given bending strength, and a given thickness ratio. For each of the 
structural conditions the case of a thin— skin structure and a solid- 
section structure is considered, since the optimum airfoil profile may be 
expected to depend somewhat on the type of structure. Attention is 
focused on the fact that the basic idea employed in the analysis involves 
the minimizing of pressure drag for a given structural requirement; the 
results obtained with this method of approach are the same as would be 
obtained if the structural characteristic were maximized for a given value 
of the drag . 4 

Mathematical Formulation of Problem 

The pressure drag c<i of an airfoil with a thick trailing edge is 
the sum of the base drag and the pressure drag of the surface forward of 
the trailing edge. Letting P be the surface pressure coefficient, y(x) 
the function defining the surface, and P^ the base pressure coefficient, 
then C(j may be expressed as 

c d'=!^ c p y dx-p b £ CD 

The problem is to determine the particular function- y(x) and the corre- 
sponding value of the trailing— edge thickness h w^iich minimizes this 

expression for a given auxiliary condition. , 

4 This statement, which appears evident from physical considerations, is - 
equivalent to Mayer* s reciprocity theorem for isoperimetric problems in 
the calculus of variations. 



NACA TN 2264 


7 


Before expressing the various auxiliary conditions in analytical 
form, it should he noted that the surface pressure coefficient P is 
regarded as a known function of the variable y* and the two parame- 
ters y ’ 0 (surface slope at leading edge), and Moo. .The actual func- 
tional form of PCySy’o'* 1 ®-) will depend on whether linearized, second- 
order, hypersonic, or shock-expansion theory is employed in calc ula ting 
surface pressures . For exampl e, if l inear theory were employed, the 

explicit expression P = 2y' / would he used; hut, if shock- 
expansion theory were employed, a more complex implicit expression involv- 
ing y ’ 0 as well as y* and Moo would have to he used. In order to 
allow various theories to he employed, the particular functional 'form of 
P(y* >y* 0 >Moo) will at present he unspecified. The equations which result 
can be applied to any of the various theories hy substituting the appro- 
priate function for P. 


Turning now to the consideration of auxiliary conditions, it is 
clear that some integral expression will he involved, since the function 
y(x) is not known beforehand. If, for example, the airfoil is a solid- 
section structure and the moment of inertia is prescribed, then the par- 
ticular auxiliary condition which y(x) must satisfy in addition to 

minimizing c d is that the integral f y 3 dx be constant. A differ- 


ent auxiliary condition would, of course, be represented by a different 
integral. In the present investigation a somewhat generalized auxiliary 
condition is used which is represented by the single integral 


I = — / — 87 - rt dx = - constant ( 2 ) 

0 J o (t/2)° 

where n and 0 are constants. Thus the example just cited is a spe— * 
cial case of the above integral with n = 3 anfl o=0. To illustrate 
further, the auxiliary condition of a given section modulus of a solid- 
section airfoil is represented by the special case n = 3 and 0 = 1 . 

The corresponding solution for y(x) in this latter case would provide 
the profile of least pressure drag for a given bending strength. 

Some of the different structural criteria to which the general inte- 
gral ( 2 ) corresponds are summarized in the following table : 5 


5 For thin— skin structures the thickness of skin is taken to be c onstant 
over the chord length. The two cases n = 0 = 0 and n = 0 = 1 are 
not included in this table as they apparently represent no sensible 
practical problem. • 
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n 

a 

Structural criteria 

1 

0 

given torsional stiffness, or torsional strength, or 
volume of thin— skin structure 

2 

0 

given bending stiffness of thin— skin structure 

3 

0 

given bending stiffness or given torsional stiffness 
of solid— section structures 

2 

1 

given bending strength of thin— skin structure 

3 

1 

given bending .strength of solid— section structure 


Thus, by solving the problem with the general integral (2) left in 
terms of n and o a wide variety of auxiliary conditions can be obtained 
simply by substituting appropriate integers for n and or. From an engi- 
neering viewpoint the general form of equation (2) enables approximate 
solutions to be obtained for wings of intermediate structural, solidity by 
interpolating between the solution for essentially zero solidity (thin— 
skin structures) and the solution for complete solidity ( solid— section 
structure). 

Summarizing, the problem formulated can be stated mathematically as 
that of finding the airfoil— ordinate function y(x), and the trailing— 
edge thickness, h = y(c), which minimizes the drag expression (l) for a 
given constant value of the structural integral (2). The boundary con- 
ditions imposed axe that y(0)=0 . and that is given. If t did not 

appear in equation (2), this mathematical problem would be a relatively 
simple isoperimetric problem in the calculus of variations. The occur- 
rence of t, the maximum value of y(x), complicates matters because it 
is not known beforehand and, in fact, is one of the quantities to be 
determined from the given values of M ro , P^, and I. Actually, all equa- 
tions necessary for solving the problem formulated could be obtained 
directly from advanced treatises on the calculus of variations since the 
problem is a special case of the so-called "problem of Bolza with varia- 
ble end points." (Such a procedure would lead directly to equations (ll), 
(l2), and (13)). However, the necessary equations can also be obtained 
by the simple methods employed here. 

Method of Solution 


Since the pressure drag of the optimum airfoil, by definition, is 
the least possible of all airfoils having a given value of the structural 
integral (2), it fqllows that the pressure drag of any "varied" airfoil. 


6 As a first approximation the torsional stiffness of a thin solid— section 
profile is taken to be .proportional to the moment of inertia about the 
chord plane. 
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having ordinates and- slopes everywhere close to those of the opti mum air- 
foil, must he in the neighborhood of a minimum. Hence, by considering 
only infinitesimal changes 5y in the ordinate ^of the optimum profile, 
the corresponding increment in drag 8c ^ of such a varied. profile can 
be equated to zero. Since y(x) is to provide the true minimum, the 
resulting equation must hold for an arbitrary ordinate change 8y vary- 
ing with x, or for an arbitrary change in airfoil thickness 5t, or 
for an arbitrary change in trailing-edge thickness 6h, or for any com- 
bination of variations thereof, provided only that the integral (2) is 
constant for all such variations. Application of this basic principle, 
as will be seen, leads to a sufficient number of equations to determine 
the complete geometry of the optimum profile. 

A sketch of the type of optimum profile to be analyzed and the cor- 
responding varied profile is shown in figure 2. Various quantities which 
appear often in the subsequent analysis are illustrated in this figure. 

It is to be noted that allowance is made for the possibility that the 
optimum profile may have a straight midsection of length l, the opti- 
mum value of which must be determined from the analysis. The varied pro- 
file is selected such that it does not change the ordinate or the surface 
slope at the leading edge. Introducing the definition f = Py* for the 
purpose of brevity, and equating the drag of the optimum profile to the 
drag of the infinitesimally varied profile, yields 


i .<* - i' c ( f+ l^ Sy ') to_Pb ^ (3) 

' The small change in slope 8y’ is equal to — (6y), so equation (3) 
can be written ^ • 



dy' dx 



( 4 ) 


Integrating by parts. 


0 
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The chordwise distribution of the perturbation by is not entirely 
arbitrary; it must be such that the auxiliary condition is satisfied, 
namely, the value of I for the optimum profile must be equal to that of 
the varied profile. 



! = 1 r c (y + 6y) n 
c Jo (t/2+8t/2) G 


( 6 ) 


Retaining only first order differences, this expression simplifies to 


or 



y 11 dx + 


F C ny*- 1 
J o 


bydx - a 


at 

TT 



dx 


(7) 


0 



by dx — 0 c 


I 



( 8 ) 


This equation must be satisfied if both terms on the right side are multi- 
plied by an arbitrary constant X. Moreover, equation (5) must be satis- 
fied simultaneously. The arbitrary character of X enables the two equa- 
tions (5) and (8) to be combined into a single equation which must hold 
for arbitrary variations in by, 8h, and 6t. 



Thus, since X is arbitrary, this single equation implies that both 
equations (5) and (8) are satisfied. 


If at this point the variation 5t is considered to be arbitrary, 
equation (9) as written would incorrectly suggest that XocI must be. zero 
Hence it is to be expected that an additional term containing 5t exists 
in the integral expression of equation (9). 'Such a term arises from the 
contribution of the straight midsection to this integral., since over this 


region, by = bt/2. Also, y= t/2 and 


— (£-') = o 

dx \oy*/ 


over this region, 


so that separation into three different ranges of integration gives the 
final equation which must be satisfied for arbitrary variations in profile 
shape : 
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The variations Sy , 5h , and. 8t can now be conducted entirely independ- 
ent of each other. Each of the bracketed terms in equation (10) must be 
zero, if the individual variations are not zero. Re-employing the defi- 
nition f= Py' , the following equations are obtained: 


for 8y j=- 0 
for 6h ^ 0 

for X5t / 0 



( 11 ) 

( 12 ) 

(13) 


The differential equation (ll), of course, results from equating to zero 
each of the two integrals in equation (10). This differential equation-, 
therefore, need be satisfied only in the two chordwise regions covered by 
the limits of these integrals, namely, in the region from x = 0 to x = s, 
and in the region from x=s+Z to x = c. (See fig. 2.) If the optimum 
airfoil has a finite length of straight midsection (e.g., AB in fig. 2), 
the differential equation (ll) need not be satisfied in this intermediate 
region. 

Fortunately, one integration of equation (ll) can immediately be made, 
thereby lowering the order of the basic differential equation to be solved. 
Multiplying equation (ll) by y* gives 





+ Iny 0-1 y« 


.= 2 y‘y«‘ 


^- + y' s y" 

By« 


<3 2 P 

by ' 2 


* x 7- (y“) 

dx 
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From this it is seen that a first integral of the basic differential equa- 
tion (ll) is ' 

y’ 2 + X = constant (lla) 

dy» 


At the point, or points, where y* = 0 the ordinate is equal to t/2. 
Evaluating the constant of equation (lla) from this consideration yields 




(lib) 


This equation, together with equation (12), equation (13), the given value 
of I, and the boundary condition y(0)=0, determines the complete geometry 
of the optimum profile. 

Method of solution for the case of a given thickness ratio .— Atten- 
tion is called to the fact that special precautions must be taken in 
applying the foregoing analysis to the auxiliary condition of a given 
thickness ratio. For this particular case 5t is zero, thereby causing 
the last term in equation (10) to vanish automatically without requiring 
equation (13) to be satisfied; equation (13), therefore, does not neces- 
sarily apply when the thickness ratio is prescribed. Moreover, equation 
(lib) also does not necessarily apply since it was assumed in the process 
of obtaining this latter equation that the optimum airfoil had at least 
one point where y* = 0. Such is not the case for the auxiliary condition 
of a given thiclaiess ratio, and hence more detailed consideration is 
required. 

The appropriate differential equation to be. used when t/c is given 
maybe obtained from equation (lla) 7 by setting X = 0. There results 


y’ 2 = constant (l4) 

ay* 

which is satisfied by any straight surface y* =. constant, regardless of 
whether linearized, second-order, or shock-expansion theory is used for P. 
The appropriate condition which must be satisfied at the trailing edge is, 
from equation (10), 

( P + y’ 1^* " P b^ 5h '= 0 (15) 

7 If equation (ll) is used there results P+y*5p/6y« = constant, which also 
is satisfied by any straight surface. The constant in this latter equa- 
tion, however, does not have the same value for both straight segments 
comprising the profile; whereas, the constant in equation (14) is the 
same for both segments. (See appendix.) 
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Here the inequality is included since Bh for the case of a given thick- 
ness ratio is not always entirely arbitrary. Thus, when h = t/2 (wedge 
airfoil) 8h is restricted to always he negative, and a minimum c an 

/ • \ 

exist if Fi+ yi*^ by* J ~ ^ ^ 0; this would make Bc^ always positive 

instead of just making stationary. Consequently, under certain con- 

ditions two solutions are possible. First, the upper half of the optimum 
profile may consist of two straight segments with h <t (as illustrated 
in fig. 3), provided the equal sign-in (15) applies. Second, the optimum 
profile may be a wedge profile with h = t, provided the inequality sign in 
(15) applies. If both types of solution are physically possible in a given 
case, the true solution, of course, would be the one with lower drag. 

If it is possible to obtain a general solution to equations (lib), 

(12), and ( 13 ) (such is the case for linearized supersonic flow), then 
the first of the two solutions mentioned above may be obtained without 
solving equations (14) and ( 15 ), but by passing the general solution to 
the limit as n <». 8 In order to verify that this limiting condition 
represents the auxiliary condition of a given thickness ratio, three 
conditions must be satisfied: First, as n — > » the auxiliary 

integral must correspond to the case of a given thickness ratio; second, 
the differential equation (lib) must reduce to equation (l4) in the 
limit as n — > »; and third, the infinite value of n must be com— 
patable with equation (13). Consider for the time being that the chord 
is of unit length. Since for any reasonable airfoil y< 1, it follows 


y 11 dx 


and hence 


A solution for 


1=0 would represent the optimum airfoil determined without regard for 
an auxiliary integral. Such is the condition that would be used in 
determining the optimum airfoil for a given thickness ratio; hence the ' 
first of the above-mentioned conditions is satisfied. Inasmuch as 
y/(t/2)<l, it is evident that [y/(t/2)] n — > 0 as n — > », thus 
reducing the differential equation (lib) to the form 




= constant 


which is the correct differential equation. As regards the third con- 
dition to be satisfied, it is seen that the values I = 0, n = eo, 1=0, 
and a = finite are compatable with equation ( 13 ) - Consequently, the 
limiting case n — > 00 in the general solution to equations (lib), 

( 12 ), and ( 13 ) represents one of the possible solutions for the case of 
a given thickness ratio. This fact will be used later in the report. 


^t should be noted that the value n = 0, when substituted in the 
auxiliary integral, gives I = (t/2) -0 = constant; but this value 
cannot be used to obtain the solution for the case of a given thick— 
ness ratio because n = 0 1 s incompatible with equation ( 13 ) . 
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Qualitative Results Obtainable Without Specializing to 
a Given Type of Supersonic Flow 


Although few quantitative results can be obtained from the basic 
system of equations (ll), (12), and (13), without specifying a particular 
form for the surface pressure coefficient P, there is’ one general quali- 
tative result that can be obtained from equation (13) without any further 
calculation. The optimum length of straight midsection l always is 
zero for the auxiliary condition of a given torsional stiffness or a 
given bending stiffness (a=0), but never is zero for the auxiliary condi- 
tion of a given bending strength (a=l). Since (t/2) n-a would be the 
value of I for a rectangular— bar airfoil of the same thickness as the 
optimum airfoil, it is evident that for an actual airfoil l/(t/2) n-0 
will be of the order of one— half or two— thirds. As an example, this 
means that when bending strength is critical in a thin— skin structure 
(n=2) , the optimum length of straight midsection will be of the order of 
one— third the chord length. 


CALCULATION OF OPTIMUM PROFILES USING LINEAR AIRFOIL THEORY 


Specialization and Solution of General Equations 
for Linearized Supersonic Flow 


According to the theory of linearized supersonic flow, the local 
surface pressure coefficient on an airfoil is given by 



For this approximation the basic differential equation (lib) becomes 

Y ,2 | = *. [(t/2) n -y n ] (17) 

or, after solving for dx/dy. 



This can be put into a more convenient form by introducing the d ime nsion— 

less variables X, Y, L, and H defined in the list of symbols, and 

eliminating by evaluating equation (l8) at x = s. Between the 
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leading edge and the first downstream position of maximum thickness, 
dy/dx is positive while the dimensionless variables X and Y both 
vary from 0 . to 1. Along the length of straight midsection Y = 1, 
while X varies from 1 to 1 + L. Along the downstream portion of 
curved surface, dy/dx is negative and Y varies from 1 to H, while X 
varies from 1 + L to c/s. Consequently, equation (l8) 'giving .X as a 
function of Y becomes (with the. convention that the sign of all radi- 
cals is positive) 


1 - 


i r x dY 


n JY J~l-Y n 


on surface facing upstream 


X = { 


(19) 


1 + L + / dY on surface facing downstream 

k nJy J 1— Y n 


where the constant k n depends only on n and is given by the definite 
integral 



1 2 for n = 1 
it/2 for n = 2 
1.4023... for h = 3 


( 20 ) 


It may be noted here that integrals of the type occurring in equa- 
tion (19) also occur at numerous places in the subsequent analysis. Such 
an integral, being a function of the parameter n and the lower limit Y, 
can be evaluated either graphically or analytically. An analytical eval- 
uation for the first three integer values of n yields 


dY 


J 1— Y n 


2 J±-Y for n = 1 

rt/2 — arc sin Y for n = 2 

3 -1 / 4 F (k,cp) for n = 3 


( 21 ) 


where 


v • 5* 

k = sin — 
12 


0.9659... 


-fa- 1+Y 
V3+1-Y 


cos cp. = 
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The function F(k,cp) is the incomplete elliptic integral of the first 
kind of modulus k and amplitude cp. (A table of this function is given 
in reference 7> page 122.) For convenience the various formulas devel- 
oped later are left in terms of the above integral; specialization to the 
individual functions indicated in equation (2l) for a given n could be 
made in any subsequent formula if desired.- 

The results of the analysis thus far involve the same mathematical 
functions as occur in the analyses of Puckett and Drougge. Thus, their 
results could be obtained at this point by setting cr = 0, L = 0, and 
imposing the boundary conditions y(o) =0 = y(c). In the present 
analysis, however, equations (12) and (13) which determine the optimum 
trailing-edge thickness and the optimum length of straight midsection 
must now be used. For linearized supersonic flow equation (12) becomes 



or, on using the relation -k n dX = dY/^fT-Y* 1 which applies to the sur- 
face facing downstream (equation (9)), there is obtained 


PbP = g k n J l-a n 

t/c (s/c) 


( 22 ) 


As defined earlier, H = h/t is the optimum trailing-edge thickness ' 
expressed as a fraction of the maximum thickness. 


Equation (13) can be written in terms of X 

I 


and Y as 


s 


a 

n 



Y n dX + Z 



(13a) 


or, after specializing to linearized flow. 



p 1 Y* 1 dY 

o 


+ k n L + / 


pi yn dY 


Jr «/l — Y n 


(23) 
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Equations (22) and (23) can be put into more usable forms by noting from 
equation (19) that 


£ 

s 


1 + L + — 
k n 



(24) 


After some algebraic manipulation involving integration by parts and 
introduction of the definition B = -PbP/(t/c) there results from com- 
bining equations (22), (23 ) > and (24), 


2n( n+2-a ) J 1-H n 

It , 2aH Jl-E* + 

1 

*3 

H 

(n-c)(n+2) 

n(n+2-a) J 

H Jl -IP- J 


(25) 


and 


L = 


20 

k n (n-cr)(n+2) 



+ H 



+ 


r 1 dY 

J E J i_yn 


) 


(26) 


These latter two equations are the final equations determining the opti- 
mum dimensionless trailing-edge thickness H and the optimum dimension- 
less length of straight midsection L. The corresponding equations 
involving the given value of I can be developed from equations ( 13) , 
(22), (25), and (26) as follows: 



This la’st equation determines H as a function of l/(t/2) n ~ a , or vice 
versa. 


Solution for the case of a given thickness ratio .— Since the use of 
linearized theory provides a general solution in closed form of the basic 
equations (lib), (12), and (l3)> the optimum profile for a given thickness 
ratio can be obtained, according to considerations presented earlier, . 
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simply by letting n — ^ co in the general solution. Since T is less 
than unity, it is evident that, for very large values of n, 



(25), and noting that L = 0 


r 


and k n — > 1. Using equations (19) and 
for the present case, it follows that 


X = 


T on surface facing upstream 

< 

2— Y on surface facing downstream 


(28) 


2(2-H) 

(29) 

l/(2-H) = 2/B 

(30) 


Equation (28) shows that the optimum surface has a discontinuity in slope 
at X = 1, and that both segments make a common angle with the chord 
plane. Equation (29) provides the required relation between the base 
pressure parameter B and the optimum trailing-edge thickness. Equa- 
tion (30) determines the position of maximum thickness. As H varies 
from 0 to 1 equation (29) covers only the range of B from 2 to 4. 

Within this range the above equations apply, and the optimum profile is 
of the type illustrated in figure 3. For the range of B from 0 to 2 
the second possible solution discussed earlier, namely, a wedge profile, 
represents the optimum section. For values of B greater t han 4 the 
double-wedge airfoil with a sharp trailing edge is the optimum for a given 
thickness ratio. 

It is remarked that the above solution also can be obtained quite 
easily by solving equations (l4) and (15) directly, instead of employing 
the limiting process. This. direct method can be used to determine the 
optimum profile in those cases where the general solution in terms of n 
and (T cannot readily be found. Such is the case when shock— expansion 
and similar higher-order theories are employed. 


Results for Linearized Flow and Discussion 


Significance and physical meaning of the base pressure parameter .— 

A relationship that greatly simplifies the determination of an optimum 
profile in linearized flow can be deduced from the structure of equa- 
tions (25), (26), and (27). '.For a given structural requirement the dimen- 
sionless optimum profile, ' Y(x) depends only on the base pressure param- 
eter B = — ! Pj 3 P/(t/c), and not on the individual values of P^, Moo, or t/c. 
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This fact may he seen by noting from equations (26) and (27) that both L 
and l/(t/2) n-a depend only on H> vhich, in turn, depends only on B 
through equation (25). Thus, although the Mach number, base pressure, 
and airfoil thickness ratio each indirectly affect, the optimum airfoil 
profile, it is only necessary to know the value of B, the base pressure 
parameter, in order to determine the dimensionless thickness distribu- 
tion (2y/t) of the optimum airfoil section. Knowledge of both I and 
B, of course, is sufficient to determine t/c as well as the dimension- 
less distribution of thickness. 

A simple physical interpretation of the base pressure parameter can 
be given if it is recalled that the basic means by which a thickened 
trailing edge reduces the over-all pressure drag is through a decrease 
in pressure foredrag at the expense of a smaller increase in base drag. 
Thus the optimum dimensionless, distribution of thickness must depend 
essentially on the ratio of base drag to pressure foredrag. The base 
drag for a given H is proportional to (— P^)(t/c); whereas the pressure 
foredrag for a given T distr ibution is, according to linearized theory, 
proportional to ( t/c) 2 / <v /M 0( f— 1. 

Hence, 

. base drag -Pb(t/c) = -Pb,/M</-1 

pressure foredrag ~ (t/c) 2/yjjQ ° t /c ' “ 


or, in words, the base pressure parameter is proportional to the ratio 
of base drag to pressure foredrag. 


Condition under which optimum profile has a-blunt trailing edge .— 
From equation (25) it is easy to deduce the condition under which the 
optimum airfoil will have a blunt trailing edge. The critical condition 
is obtained by setting H = 0. This determines a particular value of B, 
say B . 


B 


r 


4nk n ( n+2-0 ) 
(n-a)(n+2) 


8 

6.283.. . 

5.609.. . 

9.425.. . 

6.730.. . 

4 


for n. = 1, 0 = 0 

for n = 2, a = 0 

for n = 3> cr = 0 

for n = 2, 0 = 1 

for n = 3) a = 1 

for n = oo, 0 finite 


(31) 


A lower value of B would correspond, for example, to a lower base drag, 
hence the physical significance of B can be stated' quite simply: the, 
optimum airfoil has a blunt trailing edge for B <B; whereas it has a 
sharp trailing edge for B >B. Excluding the £ase of a given thickness 
ratio (n = 00, 0 finite), the average value of B for the various n and 
0 itemized in the above equation is a little greater' than 7. This 
leads to the following approximate result: If the base pressure parameter 
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B s -PbP/(t/c ) is less than about 7> then the optimum airfoil for a given 
structural requirement has a blunt trailing edge. 

An analogous result that may be of some interest concerns the opti- 
mum profile for the extreme condition when a vacuum exists at the base. 

By setting = Pjj V =— 2/7M 0 f in the equation defining B, and using 
equation (31), a relation is determined between Mach number and airfoil 
thickness ratio. Employing’ a bar to denote the critical condition for 
H = 0, and a subscript v to denote a vacuum at the base, there is 
obtained 

c) v 

With the approximation B as 7 there results for the limiting thickness 
ratio beyond which the optimum profile has a blunt trailing edge even 
though a vacuum exists at the base: 

0.10 for M,*, = 1.5 

.09 for M«, = 2 

(t/c) v ~ / .06 for Mco = 3 (33) 

.05 for Moo = 4 

.02 for Moo = 10 

With a thickness ratio of 5 percent, for example, it is evident that 
regardless of the base pressure the optimum airfoil has a blunt trailing 
edge at any Mach number beyond about 4. From the above numerical values 
it also may be inferred that the optimum trailing-edge thickness increases 
if either the Mach number or the airfoil thickness ratio is increased. 

This fact is more apparent from various examples presented later. 

Comparison with results of other investigations .— As a partial 
check on the equations developed, several limiting cases can be obtained 
by specializing to particular values of n, a, and H. First, if the 
base pressure coefficient is zero, corresponding either to zero base drag 
or else infinite Mach number, then B = 0. From equation (25) it follows 
that H = 1. In other, words, the optimum profile for F^ = 0 has its 
maximum thickness at the trailing edge. If the Mach number is finite 
and the base drag zero, then this result checks simple physical consider- * 
ations. If the Mach number is infinite (for which B = 0 even if a 
vacuum exists at the base), then this result checks the qualitative core- 
sideration of Saenger referred to in the introduction. 

A second limiting case that easily can be checked may be obtained 
by considering only the auxiliary conditions of given stiffness of sharp- 
trailing— edge profiles. The appropriate results are obtained by setting 
0=0 anid H = 0. From equation (26) it follows that 1 = 0 . From 
equations (19) and (2l) it is seen that the optimum sharp— trailing— edge 
profile is a doubly symmetric profile, each side of which is comprised of 
the arc of a parabola for n = 1, the arc of a trigonometric sine function 


(t/ 
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for n = 2, and the arc of an elliptic sine function for n = 3; these 
are the results obtained previously by Puckett ana Drougge. 

Summary curves of the principal results .— In figure 4 the optimum 
dimensionless trailing-edge thickness H is plotted as a function of the 
base pressure parameter B. Each curve in this figure is obtained by 
substituting the indicated values of n and a in equation (25). It is 
to be remembered that the curve consisting of three straight-line segments, 
corresponding to n = » and o = finite, represents the auxiliary condi- 
tion of a given thickness ratio. The other values of. n anri a repre- 
sent the various structural criteria listed in the table, page 8. 

The location of the optimum position of maximum thickness s/c, as 
determined by equations (24) and (26), is plotted in figure 5 as a func- 
tion of B. The values of n and a used here are the same as in 
figure 4. Comparing these two figures it can be seen that, as would be 
expected, the optimum position of maximum thickness moves steadily rear — 
ward as the optimum trailing-edge thickness is increased. 

Curves relating the value of I to the base pressure par ame ter are 
shown in figure 6. These curves represent equation (27). Since I is 
related to the optimum length of straight midsection through equation (13), 
the ordinate in this figure represents either of the two equal quantities, 
l/(t/2) n-cr or nZ/ac. Figure 6, therefore, can also be used to determine 
Z/c in those cases where a is not zero. If a is zero, then Z is 
zero, as noted before. 

The curves in figures 4, 5, and 6, together with the fact that all 
segments (except the straight midsection) of the optimum profile are 
composed of certain characteristic curves, provide sufficient data for a 
. rapid determination of the optimum profile. These characteristic curves 
previously were shown to be arcs of a parabola, a trigonometric sine 
curve, an elliptic sine function, or a straight line, depending respect- 
ively on whether n is 1, 2, 3, or » . 

Graphical method of determining an optimum profile .- The dependence 
of all dimensionless geometric quantities, such as H, s/c, and Z/c, 
on the single dimensionless parameter B enables a simple and concise 
graphical method to be formulated for determining an optimum profile. 

This method is illustrated in figure 7(a) for four auxiliary conditions. 
Attached to each curve in this figure is a scale unevenly divided in 
units of B. Since B determines H (fig. 4), the uneven divisions on 
each B scale are so' located that. a perpendicular dropped from each divi- 
sion intersects the curved arc at the point where the local ordinate is 
H times the m axim u m ordinate. The graphical procedure of determining 
an optimum profile consists of one step; namely, of drawing a line which 
is perpendicular to the chord, and which, if extended, would pass through 
the appropriate point on the B scale. Of the two areas so defined, the 
one which is cross sectioned in figure 7(a) represents the upper half of 
the optimum profile. To illustrate this procedure for the auxiliary 
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conditions considered in figure 7(a), the example employed has arbitrarily 
been selected to correspond in each case to B = 3. In order to obtain 
the true geometric profile, ordinates of the cross— sectioned area would 
have to be multiplied by a constant factor. This factor, of course, must • 
correspond to the particular thickness ratio used in deter min i n g B. 

As is evident from figure 7(a), each auxiliary condition considered 
therein corresponds to a/n =0. If a/n is not zero, then the graphical 
procedure for determining an optimum profile is slightly more complicated. 
This is illustrated in figure 7(b), where in each case two B scales are 
employed; one for determining the optimum trailing-edge thickness, as 
before, and a new one for determining the optimum length of straight mid- 
section. The new B scale is determined from the curves of figure 6 in 
the same manner as the original B scale is determined from the curves 
of figure 4. Thus, two lines perpendicular to the chord are drawn such 
that their extensions pass through the appropriate positions on the B 
scales. The upper half of the optimum profile is obtained by joining the 
two cross— sectioned areas so determined. 

Some examples of optimum profiles .— Using the graphical method out- 
lined above, several examples have been determined for the case n = 1, 
a = 0, in order to visually illustrate the effect of Mach number and 
airfoil— thickness ratio on the optimum profile. The effect of Mach num- 
ber is shown in figure 8, and the effect of thickness ratio in figure 9. 

In both figures the base pressure coefficient has been taken as 0.6 of 
that for a vacuum. This value corresponds roughly to the experimental 
base— pressure values for a wedge airfoil reported in reference 8, and 
also to the value found at the highest Reynolds number investigated in 
reference 9* The use of the value Pb/Pbv =0.6 is sufficiently real- 
istic for purposes of illustration, although more accurate values would 
have to be used for engineering purposes. A lower value of Pb/^bv 
would result in a greater trailing— edge thickness than that shown in the 
examples, and, conversely, a higher value would result in a smaller 
thickness. 

It is apparent from figures 8 and 9 that above a Mach number of 
about 1.4 the effect of an increase in Mach number is the same as an 
increase in airfoil— thickness ratio - both result in an increase in the 
optimum dimensionless trailing-edge thickness. Below a Mach n umb er of 
about 1.4, however, the linearized theory indicates the opposite trend 
as far as the effect of Mach number is concerned.- This is illustrated by 
the top two examples and the bottom example in figure 8, and is a result 

which would be expected since the parameter B = (— P^) /M^-l/(t/c) van- 
ishes at M,^ 1 as well as at ‘M^oo. The existence of this trend in an 
actual flow at low supersonic Mach numbers is dubious because of the inap- 
plicability of linear airfoil theory near Mach number of one. 

The trends illustrated in figures 8 and 9 can be explained from, 
elementary physical considerations if it is recalled that B corresponds 
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to the ratio of base drag to pressure foredrag. Thus, H approaches 
unity as M approaches unity because the pressure foredrag in linear 
theory approaches infinity while the base drag remains finite. Moreover, 
H also approaches unity as M approaches infinity because the base drag, 
proportional to l/M 2 becomes small compared to the pressure foredrag, 
which in linear theory becomes proportional to l/M. By the same token, 

H approaches unity for very thick airfoils because the base drag, pro- 
portional to t/c, again becomes small compared to the pressure foredrag, 
which in linear theory is proportional to (t/c) 2 . 


Drag coefficient of optimum profiles .— Since the dimensionless 
thickness distribution of an optimum profile is completely determined by 
the base pressure parameter, it is to be expected that the q uant ity 
3c d /(t/c) 2 also will depend only on B. From equation (l) and the defi- 
nition of B it is seen that 


Pea 

(t/c) 2 


23 c 


P y’dx + B ; 

"C 


Substituting P = 2y‘/f3 and changing to the dimensionless variables X, 
X, and H yields 


(vfc ■ « ♦ bh 

The integral can be expressed in terms of l/(t/2) n-a . by noting that 

C /g 

(dX/dX) 2 = kn^l-Y 11 ) and that l/(t/2) n ^ = (s/c )/ Y n dX. There 

Jo 

results for the pressure drag coefficient of the optimum profile in 
linearized flow: 


Pcd = .2 1- [l/(t/2) n ~ g ] 

(t/c) 2 (s/c) 2 


+ BH 


(34) 


Inasmuch as H, s/c, ^and l/(t/2) depend only on the base pressure 
parameter, the quantity on the left side of equation (34) also depends 
only on B for given values of n and cr # Xn an actual case the drag 
coefficient can be readily calculated from equation (3^-) by using figures 
4, 5, and 6 to determine H, s/c, and l/(t/2) n-a as a function of B. 


It is of interest to compare the pressure drag of the theoretically 
optimum profile with that of more conventional sharp— trailing— edge pro- 
files. According to linear theory the drag coefficient of a biconvex 
circular-arc airfoil (c dca ) of thickness t ca is given by 
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P( c d) ca = 16 

(tea/ 0 ) 3 


A calculation of the value of the auxiliary integral for a circular-arc 
profile (l C a) I s readily made by substituting Y=2(t ca /c)x [1— (x/c)] in 
equation (2). It is found that 


"ca 


2 511 (nlT 
(2h+l)l 


(^ca/^) 


n-0 


By requiring that I ca = I where I is the value of the auxiliary 
integral for the optimum profile. of thickness t and position of maxi- 
mum thickness at s/c, then equations (34) and (35) can be divided to 
yield 



_l/(t/2 )** 

| +BH (s/c) 2 


f (s/c) 2 

{ 

' (2n+l)j 

? ] [ i/ft/s)”^ 

1 — 

^ n-a 


(36) 


This equation gives the ratio of the pressure drag of an optimum profile 
to that of a sharp— trailing— edge, circular-arc profile having an equal 
value for I. If the pressure drag coefficient of a double-wedge profile 
( c d)dw ■*’ s use( ^ as a basis of comparison instead of a circular-arc pro- 
file, there results in a similar manner 


c d 

^ c d^dw 


n 


1 - 


l/ ( t/ 2 ) 


n-cr 


| + BH (s/c y 


4(s/c) 2 j"(n+l) 


l/(t/2) 


re-cr 


2 

n— cr 


(37) 


It may be noted that the right side of equations (36) and (37) depend 
only on the base pressure parameter if the values of n and a are 
given. 


In figure 10 the two quantities c^/(c d ) ca and c^/ ( c^)^ w 
are plotted as a function of B for various values of n and a. 
Depending bn the value of B, it is apparent that the pressure drag of 
the optimum profile may be anywhere from a few percent to as much as 
75 percent less than the drag of an equivalent circular-arc sharp— 
trailing-edge airfoil. The structural criterion for which the greatest 
drag difference exists is that of a given bending strength of a thinr- 
skin structure (n=2, a=l). As an example of conditions under which 
the pressure— drag reduction is unusually large, the following conditions 
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may be cited: Mco = 2, t/c = 0.075 and Pb = -0.08. (This value of Pb 
is the experimental value reported in reference 9 for R = 0.2 X 10 s .) 
These particular numbers correspond to B = 1.85, and for the case of 
n = 2, a = 1, the relative pressure drags are c< 3 ./(c<i) ca = 0.40 and 
c d/ ( c d^dw = ® .21. Other. conditions, besides low Reynolds numbers, under 
which large pressur e-drag reductions can be obtained by employing an 
optimum blunt— trailing— edge profile are for flow at -high Mach numbers, or 
for the flow about relatively thick airfoils. The curves of figure 10(b) 
clearly illustrate the high drag of a double-wedge profile when it is 
compared to the optimum profile on the basis of a given structural 
requirement. 

If it is desired to compare the optimum profile with a corresponding 
sharp-trailing-edge profile on the basis of relative I for a given c d , 
rather than on the basis of relative ca for a given I, then the fore- 
going calculations can be applied by making only minor modifications. 

As noted earlier, the thickness distribution of the optimum profile having 
maximum I for a. given ca is the same as that of the optimum profile 
having a minimum ca for a given M. By using the subscript s to denote 
a sharp— trailing— edge airfoil (e .g. , biconvex, or double-wedge), and no 
subscript to denote the optimum profile, the relation 

n-o 

: d =const. ' a ' / I=const. 

can be deduced if it is remembered that I varies as the (n — a) power 
of the thickness, and that the pressure drag in linearized theory varies 
as the square of the thickness. The above equation shows that in employ- 
ing an optimum section the relative structural improvement that can be 
obtained for a given drag is related in a simple way to the relative drag 
reduction 'that can be obtained for a given structural requirement. 



CONCLUDING REMARKS 


The general method presented for computing the profile shape hav ing 
minimum pressure drag at zero lift has been developed for the auxiliary 

condition that I a | f [y I1 /(t/2) a ] dx is constant. For a given air- 

foil theory, the determination of an optimum profile under this condition 
involves the simultaneous, solution of equations (ll), (12), and (13), 
which are general in that the surface pressure coefficient P(y« ) and 
the parameters n and a are arbitrary. Such generality is useful since 
it allows either linear theory, second— order theory, or shock— expansion 
theory to-be used in determining the optimum profile for a number of 
practical auxiliary conditions, such as prescribed bending strength or 
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given torsional stiffness. As an illustration of the method, a solution 
has been developed in detail using linearized flow, that is, using the 

expression P = 2y , / v /M 00 2 - 1. In this simple case a complete solution in 
closed form is obtained for the thickness distribution of the Optimum 
profile . 

The principal result of the analysis for linear supersonic flow is 
that the dimensionless thickness distribution of the optimum profile 
depends only on the single parameter B= — Pfc a/M^-I/ ( t/c ) . This 
parameter has been termed the base pressure parameter, and has a simple 
physical significance in that it is proportional to the ratio of base 
drag to pressure foredrag. The dependence of an optimum profile in 
linear flow on one parameter only enables summary curves to be plotted 
showing all principal results as a function of B (figs. 4, 5, and 6). 
These curves form the basis of the graphical method formulated for 
determining an optimum profile (fig. 7) • Other results of interest con- 
cern the individual effect on the optimum profile of varying the base 
pressure coefficient, airfoil— thickness' ratio, and Mach number. The 
optimum dimensionless trailing-edge thickness increases if either the 
base pressure is increased, the airfoil— thickness ratio is increased, or 
the Mach number is increased beyond about 1.4. 

At low supersonic Mach numbers the theoretical results obtained are 
questionable since the assumptions of linearized airfoil theory break 
down as the Mach number approaches unity. The results can be applied 
safely only to cases where linear theory satisfactorily predicts the 
pressure foredrag. Although at high supersonic Mach numbers the results 
obtained under the. assumption of linearized flow also would not be. 
expected a priori to be of quantitative value, they predict, nevertheless, 
the correct result that the optimum trailing-edge thickness for infinite 
Mach number is equal to the maximum airfoil thickness. In view of this 
exact agreement in the extreme case it is conjectured that the linear 
theory fortuitously may provide a reasonable est ima te of the optimum 
trailing-edge thickness for any supersonic Mach number not close to unity. 
As regards the optimum profile shape forward of the base, however, such 
fortuitous conditions cannot be expected, since the linearized approxima- 
tion at high Mach numbers overestimates the suction forces and underesti- 
mates the positive pressure forces. This causes the calculated optimum 
profile to have too large a leading-^edge angle, a position of maximum 
thickness too far forward, and too small an inclination of the surface 
behind the position of maximum thickness, (in reference 10 some calcu- 
lations using second— order theory are presented which illustrate this 
effect on the optimum sharp-trailing— edge profile for the auxiliary con- 
dition of a given thickness ratio.) 

Because the optimum profile, by definition, has the least pressure 
drag possible under given conditions, small changes in profile shape 
would result in second— order changes. in drag. This allows some flexi- 
bility in modifying the theoretically optimum profile to more closely 
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suit individual design requirements, and means that it is not important 
to rigorously adhere to the exact parabolic, trigonometric— sine, or 
elliptic-sine contour (provided, of course, that the end points of the 
modified contour are located approximately in the optimum positions). It 
is important to adhere reasonably close to the calculated optimum trailing- 
edge— thickness, since this quantity can greatly affect the drag. In par- 
ticular, a trailing-edge thickness considerably greater than the optimum 
should not be used. Excessive trailing-edge thickness at low and moderate 
supersonic Mach numbers will result in an increase in drag rather than a 
decrease. 

The possible effects of viscosity on the theoretically optimum pro- 
file for minimum pressure drag are considered in the analysis only inso- 
far as they affect the base pressure. Shock-wave boundary— layer inter- 
action is known to modify considerably the calculated pressures near the 
trailing edge of a sharp-trailing-edge airfoil. The optimum blunt- 
trailing— edge airfoil, however, always has a lower pressure on the base 
than on the surface just upstream of the base. (For linearized flow 
equation ( 12 ) shows that Pb-Pi=Pi~Poo. ) Because of the resulting expan- 
sion, viscous effects probably are unimportant forward of the blunt trail- 
ing edge. 

It is to be remembered that some caution should be exercised in 
applying the quantitative 'results of this investigation, since the calcu- 
lations consider only pressure drag; whereas, a practical case usually 
is concerned with profile drag. The results can be directly applied if 
the skin— friction coefficient of the optimum profile is insensitive to 
small changes in profile shape (turbulent flow at high Reynolds numbers), 
or if friction is small compared to the pressure drag (thick airfoils). 

One case where difficulty possibly may occur in this regard concerns the 
optimum profile for a given thickness ratio. Here the airfoil for mini- 
mum pressure drag has a sharp corner. What effect this may have on 
boundary-layer transition is not as yet known. For most auxiliary condi- 
tions though, the optimum profiles have a smoothly curving surface with 
a continuously decreasing pressure from leading to trailing edge, and 
hence would be expected to be favorable for long runs of inminn-r flow. 

Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 

Moffett Field, Calif., Oct. 3 , 1950.. 
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APPENDIX 

NOTE REGARDING DISCONTINUOUS CHANGE IN SLOPE 
OF AN OPTIMUM PROFIIE 


In the general analyses of optimum profiles for a given thickness 

ratio it was found that the function y* — necessarily was constant 

oy' 

along each straight segment of the profile (equation (l4) ) . . No informa- 
tion was obtained, however, about the relative value of this constant 
for the two segments. The required information can readily be obta ine d 
by considering the change in drag due to a change only in position of 
maximum thickness, that is, a change in slope of both straight surfaces 
illustrated in figure 3 with no change in t or h. Using subscripts 
u and d to denote- surfaces facing upstream and downstream, respect- 
ively, it follows that . 


6c d = S(c d ) u + 5(c d ) d 


Since y' u = t/s and y* d = (t-h)/(c— s), the minimizing of c fl requires 
that 


6c d = 0 = - (^) ■ .£- 8s + 5s 

w'/u s 2 W’/cL (c— s) 2 


0 = 5s 


Thus y 


«2_ 


8p 

8y« 


y’u 


8p 

8y* 


— V T 


■y d 


u 


8pA 

Sy'^d -I 


must be continuous at the corner — a result which was used 


without proof in the general analysis. It may be noted that an alternate 
proof of this result can be obtained in an extremely easy way from the 
following known result of the calculus of variations: The Weierstrass E— 

function is continuous at the point of discontinuity on a boundary. The 


E— function in the present case is y’ 2 


8 P 

dy* * 


For auxiliary conditions other than a given thickness ratio it was 
tacitly assumed in the analysis that the optimum surface everywhere had 
a continuous slope. This assumption also requires some justification. 
It is shown in the calculus of variations that at all points of free 
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df 


variation it is necessary for y— to be continuous. For the present 

+ y*^£ 


problem this means that P + y* must be continuous at all such 

dy* N 

L. 

points. According to linear theory, P + y* = — y 1 , hence, within 

by* P 

the scope of linear theory, the surface slope' y' is continuous at all 

dp 

points of free variation. For shock-expansion theory is positive, 

dy* 

and a corner would cause a discontinuous decrease in P, y*-, and 


P + y’ ; -hence, also within the scope of shock-expansion theory the 
dy’ 

.surface slope of the optimum profile is continuous at all points of free 
variation. This justifies the assumption of continuous slope employed 
in the general analysis for auxiliary conditions other than a given 
thickness ratio. 
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Figure /.- Sketch comparing a typical sharp-traiUng-edge 
airfoil and a blunt-trailing-edge airfoil of equal area 
(equal torsional stiffness for a thin -skin structure). 



Figure 2. - Sketch of upper half of varied and optimum 

profiles . 
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Figure 3.- Sketch of upper half of optimum profile for a 

given thickness ratio. 



Figure 4. - Optimum trailing- edge thickness for linearized 

supersonic flow. 
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0 2 4 6 8 10 

Base Pressure Parameter, B= - P b fit* -/ /(t/c) 


Figure 5- Optimum position of maximum thickness for 
linearized supersonic flow. 



0 2 4 6 8 10 

Base Pressure Parameter, B-- P b /A/J -/ /(t/c) 


Figure 6. - Curves relating I and 1 to the base pressure 
parameter for linearized supersonic flow. 
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Figure 7.- Curves and scales for graphical determination 
of optimum profiles for various auxiliary conditions- 
linearized flow , illustrative examples for B = 3, vertical 
scale expanded. 
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Surface facing 
upstream 


Midsection and surface 
facing downstream 


(b) Auxiliary conditions requiring two B- scales. 


Figure 7. - Concluded. 




ct of Mach 
ed flow , vertical 
6, t/c -- 0.06. 
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Figure 9.- Examples illustrating the effect of airfoil 

thickness ratio on the optimum profile-, linearized flow, 
vertical scale expanded five times, n = l, a=0, 

Pt/Pbv-0.6, M„:2. 
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Base Pressure Parameter, B = -P b 1 /(t/c) 

(a) Drag of optimum relative to biconvex circular - arc airfoil. 



0 2 4 6 8 !0 

Base Pressure Parameter, B= -P b /(t/c) 


(b) Drag of optimum relative to double-wedge airfoil. 

Figure 10. - Drag of optimum airfoils as compared to the drag 
of two different sharp-trailing- edge airfoils in linearized 
supersonic flow. Comparison made on the basis of equal values . 
for the auxiliary integral. '^Qaca^ 
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